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Abstract 

The hypothesis that aquaporins and carbonic anhydrase (CA) are involved in the regulation of stomatal (g s ) and meso- 
phyll (g m ) conductance to C0 2 was tested in a short-term water-stress and recovery experiment in 5-year-old olive 
plants (O/ea europaea) growing outdoors. The evolution of leaf gas exchange, chlorophyll fluorescence, and plant 
water status, and a quantitative analysis of photosynthesis limitations, were followed during water stress and recov- 
ery. These variables were correlated with gene expression of the aquaporins OePIP1.1 and OePIP2.1, and stromal 
CA. At mild stress and at the beginning of the recovery period, stomatal limitations prevailed, while the decline in g m 
accounted for up to 60% of photosynthesis limitations under severe water stress. However, g m was restored to con- 
trol values shortly after rewatering, facilitating the recovery of the photosynthetic rate. CA was downregulated during 
water stress and upregulated after recovery. The use of structural equation modelling allowed us to conclude that 
both OePIP1.1 and 0ePIP2.1 expression could explain most of the variations observed for g s and g m . CA expression 
also had a small but significant effect on g m in olive under water-stress conditions. 

Keywords: AQP, CA, mesophyll conductance, stomatal conductance, OePIP, O/ea, photosynthesis limitations, water stress. 



Introduction 

Water stress is considered the main environmental factor lim- 
iting photosynthesis, plant growth, and yield worldwide, espe- 
cially in semi-arid areas, where Olea europaea is well adapted 
(Boyer, 1982; Lawlor, 1995; Flexas et al, 2004). Under 
water-stress conditions, which are related to water deple- 
tion and/or high atmospheric vapour pressure deficit (VPD), 



photosynthesis decreases through several mechanisms includ- 
ing stomata closure, reduced mesophyll conductance to C0 2 
(gm), and feedback regulation by end-product accumulation 
(Nikinmaa et al., 2013). As soil water deficit and VPD result 
in reduced stomatal conductance to C0 2 (g sc ) (reviewed by 
Lawlor & Cornic, 2002; Flexas et al., 2004) and g m (Flexas 



Abbreviations: ABA, abscisic acid; AQP, aquaporin; a.r.L, after resuming irrigation; a.w.i., after withholding irrigation; CA, carbonic anhydrase; Q, internal C0 2 con- 
centration; ETR, electron transport rate; GMT, Greenwich Mean Time; g m , mesophyll conductance to C0 2 ; g s , stomatal conductance to C0 2 ; PSII, photosystem II; 
REW, relative extractable water; RWC, relative water content of leaves; SE, standard error; VPD, vapour pressure deficit. 
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et al., 2002, 2004; Centritto et al, 2003; Galmes et al, 2007a; 
Warren, 2008a; Peeva and Cornic, 2009) in many species, sev- 
eral authors have suggested a possible co-regulation of g sc and 
g m (Centritto et al., 2003; Flexas et al., 2002, 2012). However, 
such co-regulation depends on the species and the prevailing 
conditions, such as the combination of drought with VPD 
(Perez-Martin et al, 2009) or with radiation (Flexas et al, 
2009; Galle et al, 2009). 

Under natural conditions, water stress normally develops 
gradually over periods of weeks or months, and hence it is 
possible that some acclimation occurs in addition to day-to- 
day variations in response to variable environmental condi- 
tions (Flexas et al, 2006a). Acclimation to water stress may 
comprise responses involving modification of gene expres- 
sion and plant physiology and morphology, taking place over 
days to weeks, which lead to a homeostatic compensation for 
the initial negative effects of water stress on photosynthesis 
(Lambers et al. , 2008). Nonetheless, little is known about the 
acclimation of photosynthesis to water-stress conditions in 
the short term, so studies that specifically address this issue 
are needed (Flexas et al, 2006a). The recovery phase after 
relief of stress (i.e. rainfall or irrigation) becomes another 
important part of the overall plant physiological response 
to a water-stress period. The capability for photosynthetic 
recovery from a water-stress period determines the future 
growth and survival of plants in their habitat, and depends 
on the degree and velocity of photosynthesis decline during 
water depletion (Flexas et al., 2006a, 2009). Although there 
is much information on the regulatory mechanisms of the 
response of g s to water stress and recovery (Buckley, 2005), 
little is known about the regulation of g m (Flexas et al, 
2008). The main limiting factor for photosynthesis dur- 
ing water stress or recovery can vary depending on species 
(Galmes et al, 2007a; Ennahli and Earl, 2005), the intensity 
of previous stress (Flexas et al, 2009), light and temperature 
(Galle et al, 2009), plant age (Varone et al, 2012), and the 
application of successive drought and recovery cycles (Galle 
et al, 2011). 

Therefore, besides photosynthetic biochemistry, the regula- 
tion of g s and g m is crucial to understand both processes of 
acclimation to stress and recovery after stress. Despite being 
one of the most studied physiological variables, the regulation 
of stomata is not yet fully understood (Buckley and Mott, 
20 1 3). There is, however, a large consensus on the fact that two 
main components operate to produce the adequate response 
of g s to environmental stimuli: a hydropassive loop, related to 
the hydraulic capacity of the plant, and a hydroactive loop, 
linked to both the chemical signalling and the photosynthetic 
capacity of the leaf (Buckley et al, 2003). Regarding the reg- 
ulation of g m , recent evidence suggests that anatomical traits, 
such as cell-wall thickness and chloroplast distribution, are 
among its stronger determinants (Flexas et al, 2012; Tosens 
et al, 2012; Tomas et al, 2013). Nevertheless, the most likely 
candidates proposed for the rapid regulation of g m are aqua- 
porins (AQPs) and carbonic anhydrase (CA). There are sev- 
eral pieces of evidence of AQP involvement in g m regulation. 
First, inhibiting the activity of some aquaporins by HgCl 2 , 
Terashima and Ono (2002) found a decreased g m in Vicia 



faba. Secondly, the tobacco aquaporin NtAQPl facilitates 
C0 2 membrane transport when inserted in Xenopus oocytes 
(Uehlein et al, 2003). Finally, the strongest evidence comes 
from the observed altered g m in transformed and mutants 
plants (Hanba et al, 2004; Flexas et al, 2006ft; Heckwolf 
et al, 2011; Kawase et al, 2013). 

On the other hand, as aquaporins accumulate in cells 
around stomatal cavities and in guard cells themselves (Otto 
& Kaldenhoff, 2000), they may also be involved in the regula- 
tion of g s . In fact, transgenic and mutant plants with altered 
AQPs present differences not only in g m but also in g s (Hanba 
et al, 2004; Flexas et al, 2006ft; Heckwolf et al, 2011). 
Recently, Pou et al (2013) found a strong correlation between 
g s and the gene expression of some particular AQPs during 
water stress and recovery in grapevines. Despite all these find- 
ings, the potential relationship between the gene expression 
of AQPs and the variations in g s and g m is still unclear, as 
patterns of AQP expression are complex and different AQP 
forms may induce distinct responses (Tyerman et al, 2002; 
Alexandersson et al, 2005, Galmes et al, 2007ft) 

Concerning the involvement of CA in the regulation of g m 
less evidence has been found. CA might have a role in the reg- 
ulation of g m through the establishment of the dynamic equi- 
librium between C0 2 and HCO3 (Tiwari et al. , 2005; Tholen 
and Zhu, 2011; Flexas et al, 2012). Likewise, the presence of 
isoforms of CA in all the compartments composing the dif- 
fusion way of the liquid phase of C0 2 (except on the cellular 
wall) suggests that they may be involved somehow in the reg- 
ulation of g m (Evans et al. , 2009; Terashima et al. , 201 1). The 
elimination of the stroma CA might result in a 50% reduction 
in g m (Tholen and Zhu, 2011) but, in contrast, studies with 
mutant plants with reduced CA have shown little reduction 
in g m (Price et al, 1994; Williams et al, 1996). Some authors 
have claimed that the role of CA could be species depend- 
ent and become more important in species with a high C0 2 
resistance in their cell walls and low g m , as is the case of many 
drought-adapted species (Gillon and Yakir, 2000). 

Olea europaea is an excellent plant model for studying 
physiological and molecular responses to water stress and 
recovery under realistic water-stress conditions because of its 
high reputation as drought-tolerant species (Fernandez and 
Moreno, 1999; Connor, 2005), as well as its importance in 
the Mediterranean landscape and agronomical impact (Sofo 
et al, 2009; Boughalleb and Hajlaoui, 2010). As no previous 
study has related AQP and CA gene expression with all the 
potential photosynthesis-limiting factors (g s , g m and biochem- 
istry) during water stress and recovery under field conditions, 
we studied the expression of two AQP genes characterized in 
olive, OePIPl.l and OePIPl.l, and a stromal CA. The exper- 
iment was carried out outdoors, under natural conditions of 
radiation and atmospheric demand, in contrast to most up- 
to-date studies of the relationships between AQPs and CA, 
which have been performed in plants grown in growth cabi- 
nets. The specific objectives of the present work were: (i) to 
quantify the photosynthesis limitations imposed by g s , g m and 
photosynthetic capacity during the water-stress and recov- 
ery period, and (ii) to relate these changes to changes in the 
gene expression of AQPs and CA in order to expand current 
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knowledge on the possible role of the AQP and CA expres- 
sion in the regulation of g s and g m . 



Materials and methods 

Plant material and water-stress treatments 

Eighteen 5-year-old O. europaea L. var. Manzanilla plants were 
grown outdoors in 50 1 pots (see Supplementary Fig. SI at JXB 
online) at 'La Hampa', an experimental farm of the Consejo 
Superior de Investigaciones Cientificas (37° 17' N, 6°3' W, altitude 30 
m), near Seville (south-western Spain). Pots contained sandy-loam 
soil, typical of this farm and whose hydraulic properties have been 
described previously by Palomo et al. (2002). The experiment was 
performed during the summer of 2009. The climate was typically 
Mediterranean, with a mild, wet season from October to April and 
hot and dry conditions from May to September (Fernandez et al. , 
2006). Six trees were pruned by one-third of their leaf area to reduce 
leaf water demand and to produce a slower evolution of water stress. 
Branches were cut at dawn after being sprayed with water and were 
covered with plastic bags to minimize embolism risk. Pruning cuts 
were covered with wound dressing. Thus, two different water-stress 
intensities were imposed by withholding irrigation in six non-pruned 
trees (S) and six pruned trees (SP) on 26 July. Another group of six 
olive trees was maintained as well watered and was used as control 
treatment (C). Finally, after 13 d of total water withholding, S and 
SP plants were rewatered daily for the following 6 weeks and allowed 
to recover. All the physiological measurements detailed below were 
taken in leaves of a similar age, including those carried out during 
the recovery experiment. This ruled out the possibility that new, 
recently grown leaves were measured. 



Climatic variables and soil water status 

Meteorological variables were measured with an automatic weather 
station (Campbell Scientific Ltd, Shepshed, UK) located 50 m away 
from the experimental trees. Average values taken over 30min of 
net radiation, global radiation, photosynthetically active radia- 
tion, wind speed, air temperature, relative humidity, and VPD were 
recorded. 

The volumetric soil water content (0, mm) of the substrate in the 
pots was determined by the time-domain-reflectrometry technique 
using a Tektronix cable tester (Model 1502C; Beaverton, OR, USA). 
Measurements were made at 7.00 Greenwich Mean Time (GMT) 
and averaged between 0.05 and 0.20 m in depth because there were 
no differences between the values recorded at either depth within 
each treatment. 6 was transformed into soil relative extractable 
water (REW) according to Fernandez et al. (1997): 

REW = (0-6> min )/(0 max -6> min ), 

where 6 is the actual soil water content, 0 min is the minimum soil 
water content measured during the experiment (mm), and 0 max is 
the soil water content at field capacity (mm). Field capacity was 
determined to be 0.25 m 3 nr 3 and the soil water content at perma- 
nent wilting point was 0.04 m 3 nr 3 . A REW value of >1 could be 
observed for several hours after water application, especially in the 
C treatment, as during that period 6 exceeded field capacity. 



Plant water relationships 

Leaf water potential was assumed to be equal to xylem pressure 
potential at the petiole in transpiring leaves and was measured in 
six leaves per treatment at pre-dawn ( fpj) and 1 1 .00 GMT ( f) with 
a pressure chamber (Soilmoisture Equipment, Santa Barbara, CA, 
USA). The relative water content of leaves (RWC) was determined 
at 1 1 .00 GMT in opposite leaves to those sampled to measure W as: 



RWC = (fresh weight - dry weight) 

/(turgid weight - dry weight) xlOO. 

Olive leaves were weighed immediately on harvest to determine 
fresh weight. To determine the turgid weight of the samples, these 
were kept in distilled water in darkness at 4 °C to minimize respira- 
tion losses until they reached a constant weight (full turgor, typically 
after 24 h). Their dry weight was obtained after 48 h at 70 °C in an 
oven. 



Leaf gas exchange and chlorophyll fluorescence 

Leaf gas exchange in leaves similar to those used for plant water 
relationships was determined simultaneously with chlorophyll fluo- 
rescence at 1 1 .00 GMT using the open gas-exchange system Li-6400 
(Li-Cor, Lincoln, NE, USA) equipped with an integrated fluores- 
cence chamber head (Li-6400^1-0; Li-Cor). Photosynthesis was 
induced with saturating light (1800 umol nr 2 s~') and an ambient 
concentration of C0 2 of 400 umol mol -1 . 

From the fluorescence measurements, the actual photochemical 
efficiency of photosystem II (<J> Psn ) was determined according to 
Genty et al. (1989) as: 

*psii =(*"«.' -F s ')/F m , 

where F s is the steady-state fluorescence in the light [here photo- 
synthetic photon flux density (PPFD) of 1800 pmol irr 2 s~'] and F m ' 
is the maximum fluorescence obtained with a light-saturating pulse 
(~8000 umol m 2 s -1 ). As <t>p S n represents the number of electrons 
transferred per photon absorbed by PSII, the rate of electron trans- 
port (J) can be calculated as: 

J = <Dp SII xPPFD x0.5 x0.93, 

where PPFD is the PPFD incident on the leaf, 0.5 is a factor 
that assumes equal distribution of energy between the two photo- 
systems (Laisk and Loreto, 1996), and 0.93 is the leaf absorptance 
determined in five leaves per treatment using an integrating sphere 
with a portable spectroradiometer (LI- 1800; Li-Cor) and calcu- 
lating absorptance as 1 - reflectance - transmittance. The rela- 
tionship between <E>psii and the quantum efficiency of gross C0 2 
fixation (<J> C o2) was obtained by varying the light intensity under 
non-photorespiratory conditions in an atmosphere containing <1% 
0 2 (Valentini et al, 1995). The correlation between <r>psii and <I> C02 
was done under different conditions of soil water deficit and VPD 
(see Supplementary Fig. S2 at JXB online) to check for poten- 
tial errors in the estimation of g m due to non-linearity of electron 
transport rate. 

Estimation of g m was performed with the 'variable J method' of 
Harley et al. (1992): 

g m = I {C; -(T *[ / + 8(4, + Bi)] /[ J 4(4 + ty] }■ 

where A K and internal C0 2 concentration (Q) were taken from 
gas-exchange measurements at saturating light, J was estimated from 
fluorescence, the rate of mitochondrial respiration in the light (R d ) 
was assumed to be the same as measured dark respiration (Warren, 
2004), and F* was estimated at the measuring temperature by using 
the kinetic constants proposed for this parameter by Bernacchi et al. 
(2002). To estimate R d , we covered the leaves (n-6) for 30min and 
measured the net assimilation of C0 2 rates using a modified 2x3 cm 
broadleaf chamber and an integrated light source (LI-6400-02B; 
Li-Cor) at a flow rate of 250 umol air s -1 , The relationships between 
<I>psii and <I>co2' and R d were measured in three stages of the experi- 
ment: one at the beginning of the experiment, the second 7 d after 
withholding irrigation, and the third 2 d after resuming irrigation. 
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No significant differences in leaf absorptance were found through- 
out the experiment, despite changes in the RWC of leaves. 



Quantitative photosynthesis limitations analysis 

To assess the limitations imposed by water stress and recovery on 
photosynthesis, a quantitative limitation analysis of photosyn- 
thesis was conducted according to Grassi and Magnani (2005), 
with modifications. This approach requires measurements of ^4 N , 
g s , gm> and F cmax to calculate the partition of photosynthesis limi- 
tations into components related to stomatal conductance (Si), 
mesophyll conductance (MC L ), and leaf biochemical characteris- 
tics (B L ), assuming a reference treatment where maximum assimi- 
lation rate, g s , g m , and F cmax can be defined (see Supplementary 
Information in JXB Online for details of equations used). Total 
limitations (T L ) were defined as the sum of S L , MC L , and B L . As 
actual electron transport rate (ETR, i.e. fluorescence-derived J) 
is tightly coupled with V cm!ix (Galmes et al., 2007a) and should 
indeed reflect gross photosynthesis (Genty et al., 1989; Valentini 
et at, 1995), B L was calculated using ETR instead of K cmax as a 
surrogate for leaf biochemistry. There were three reasons for the 
use of ETR as a surrogate of K cmax in the limitation analysis: (i) 
under our experimental approach, it was unfeasible to build a 
number of A N -Q curves enough to estimate K cmax for every sam- 
pling date; in addition to this, when g s is low, the estimation of 
^cmax is highly prone to errors due to the low rates of leaf gas 
exchange; (ii) as actual electron transport rate is tightly coupled 
with K cmax (Galmes et al. , 2007a) and should indeed reflect gross 
photosynthesis (Genty et al, 1989; Valentini et al, 1995), ETR 
can be used instead of F cmax as a surrogate for leaf biochemistry; 
the uncertainties in the determination of F cmax reported by sev- 
eral authors (Patrick et al, 2009; Gu et al, 2010) can be avoided 
by using ETR, as has been verified and confirmed by Galle et al 
(2009, 2011); and (iii) we checked that, under our conditions and 
during the progress of water stress and recovery, there was a good 
correlation between ETR and F cmax for the range measured; this 
correlation is shown in Supplementary Fig. S3 at JXB online. In 
the current study, the maximum assimilation rate, concomitantly 
with g s , g m , and ETR, was reached under well-watered conditions, 
and therefore the C plants were used as a reference. However, as 
^4 N of the C plants increased during the experiment, presumably 
due to environmental changing conditions, the values of the C 
plants for each day were considered as the reference for the S and 
SP treatments determined on the same day. 



RNA extraction and real-time reverse transcription-PCR analysis 

The same leaves used throughout the experiment for gas exchange 
and chlorophyll fluorescence measurements were harvested at 
12.00 GMT for gene expression analysis. They were harvested at 
the same time of the day to avoid diurnal fluctuations in the tran- 
scripts (Henzler et al , 1999; Laur and Hacke, 2013). Three replicates 
per treatment and day were sampled and were frozen immediately 
in liquid nitrogen. Leaf tissues of each sample were ground to a 
fine powder with liquid nitrogen using sterile mortars and pestles. 
RNA extraction was performed on lOOmg of plant material using 
an RNeasy Plant Mini kit (Quiagen, Venlo, The Netherlands) and 
treated with RNase-Free DNase (Qiagen) to avoid genomic DNA 
contamination during RNA purification. The concentration of 
total RNA fraction in each sample was estimated by measuring the 
sample absorbance (A 260 /A 2 go) with a spectrophotometer (Lambda 
6 UV-VIS; Perkin-Elmer, Buckinghamshire, UK). RNA reverse 
transcription to cDNA was performed on 1 |ig of RNA using a 
QuantiTect Reverse Transcription kit (Qiagen). The samples were 
then stored at -80 °C until real-time PCR analysis. 

The three genes whose expression was chosen to be studied as 
influenced by water stress and recovery were: OePIPl.l (Secchi 
et al, 20076; GenBank accession no. DQ202708) and OePIP2.1 
(Secchi et al. , 20076; GenBank accession no. DQ202709), encoding 



functional water-channel proteins in olive; and the gene encoding 
an olive CA enzyme located in the chloroplast stroma (GenBank 
accession no. FN8 14304). We constructed a multiple alignment 
of the whole Arabidopsis family reported by Fabre et al (2007), 
including the olive CA. The isoform of CA of olive was con- 
firmed to be aligned with Arabidopsis stroma isoforms |3CA (see 
Supplementary Fig. S4 at JXB online) described by Fabre et al. 
(2007) to be targeted to the chloroplasts. The CA full DNA cod- 
ing sequence of Arabidopsis thaliana was retrieved from phyto- 
zome (http://www.phytozome.net) for phylogenetic analysis with 
the identified O. europaea mRNA (CA, GenBank accession no. 
FN8 14304). Alignment and phylogenetic analysis were performed 
with the seaview toolkit (Gouy et al, 2010). Multiple alignment 
of mRNA sequences was performed using the MUSCLE algo- 
rithm (Edgar, 2004). As phylogenetic inference could be biased 
by partial sequences in the alignment, as with the O. europaea 
partial mRNA in our dataset, prior to phylogenetic analysis we 
extracted the evolutionary and conserved aligned blocks from the 
alignment using the Gblock algorithm (Talavera and Castresana, 
2007). A phylogenetic tree was created by the maximum-likelihood 
approach using the PhyML algorithm (Guindon et al. , 2010), with 
the K80 (Kimura, 1980) nucleotide substitution model, 100 ran- 
dom tree starts, and an optimized ratio between nucleotide transi- 
tion and transversion (TJT V ratio). Branch support was estimated 
by an approximate likelihood ratio test based on the Shimodaira- 
Hasegawa-like procedure (Shimodaira and Hasegawa, 1999; 
Guindon etal, 2010). 

Despite the fact that both OePIPs are more expressed in roots 
and twigs, we studied them because they are also expressed in leaves 
and are actually the only PIPs that have been studied and charac- 
terized in olive (Secchi et al. , 2001b). We included the p-actin gene 
(GenBank accession no. AF545569) as a housekeeping control. 
Specific primers were designed by Applied Biosystems (Foster City, 
CA, USA). The gene expression assays consisted of a 20x mix of 
unlabelled PCR primers and TaqMan MGB probes (labelled with 
carboxyfluorescein dye). 

Real-time PCR (7300 Real-Time PCR System; Applied 
Biosystems) was performed on the samples, and relative gene 
expression was determined using the relative standard curve method 
(Applied Biosystems). As a calibrator, a sample from the experiment 
(C plant) was chosen, the calibration curves being performed with 
a 1, 10, 50, and lOOng dilution of the cDNA of the same sample, 
for each of the genes chosen to be studied plus the housekeeping 
gene. Fast PCR cycles were performed running 96-well Fast reaction 
plates with 20 ul of the reaction mix in every well and using TaqMan 
Universal PCR Master Mix (2x; Applied Biosystems). The PCR 
thermal cycling conditions used were the default ones, as indicated 
by the manufacturer's instructions. 



Statistical analysis 

Statistical analysis was carried out using Statistica for Windows 
v.6.0 software package (StatSoft). Significant differences between 
means were assessed using one-way analysis of variance applying 
Tukey's test (P<0.05). Previously, data were tested for normality and 
homoscedasticity. 

A d-separation method of path analysis was used to test several 
conceptual models that could explain the covariance among the 
measured physiological variables, especially g m , and the genetic 
expression of AQPs and CA. The proposed models were based on 
previous knowledge and the hypothesis of interactions between vari- 
ables. Causal relationships between variables were combined to form 
directed graphs (the path models). These directed graphs implied 
a series of independence relationships between pairs of variables. 
Pearson correlation coefficients between pairs of variables were 
determined. The relationships among variables were then translated 
into a structure of variances and covariances that could be tested 
against the observed data. Path analysis was performed using the 
Causal Toolbox (Shipley, 2000) packages. 
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Results 

Stress period 

REW dropped rapidly after water withholding in the S and 
SP treatments (Fig. la). Both stress treatments showed REW 
of <0.3 as soon as 2 d after withholding irrigation (a.w.i.). 
After that, a gentle decrease in REW down to 0.08 was 
recorded, suggesting very little water use by the plants. On 
average, the atmospheric demand for water was high, with 
maximum VPD values usually >4 kPa, and even reaching 6.8 
kPa (Fig. lb). During the initial drop in REW, both S and 
SP plants maintained both W and RWC values similar to the 
C plants (Fig. 2), although with slight larger values for SP. 
However, differences with C were observed on day 4 a.w.i. 
and afterwards, with higher values in SP (with a smaller leaf 
area) than in S. At the end of the stress period, both stress 
treatments reached similar minimum values of W and RWC 
of approximately -5.5 MPa and 55%, respectively. Leaf gas- 
exchange variables were more responsive to the decline in soil 
REW, and a reduction in all three variables, g s , g m , and A N , 




2 34 5 
Days a.r.i. 



Fig. 1 . (a) Soil REW (see equation in Materials and methods) throughout 
the experiment in pots of the three treatments: C, control plants; S, stress 
plants; SP, stress-pruned plants. Values on each date are the means (±SE) 
between two different depths (0.05 and 0.20 m). Different letters indicate 
significant differences among treatments within each date (analysis of 
variance, Tukey: P<0.05). (b) VPD of air during the experiment in July and 
August of 2009. S and SP plants were last irrigated on day 0 a.w.i., with 
daily irrigation being applied again 13 d later (discontinuous line), a.w.i., 
after withholdong irrigation, a.r.i., after resuming irrigation. 



was observed before changes in leaf water status were notice- 
able (Fig. 3). In parallel, ETR showed a continuous decrease 
in both S and SP plants from the beginning (Fig. 4). 

OePIPl.l expression in S and SP plants peaked on day 
4 a.w.i. and returned to C values at the end of the drought 
period (Fig. 5a). Expression of OePIP2.1 in stressed plants 
was similar to that in C plants at the beginning of the drought 
period but diminished at the end (Fig. 5b). For CA, relative 
expression in stressed treatments was similar to that in C 
plants during the first days of drought but diminished at the 
end of this period (Fig. 5c). However, compared with the first 
day, relative expression in stressed treatments showed a con- 
tinuous decrease from the beginning, similar to that described 
for ETR (Fig. 4). 



Recovery period 

After resuming irrigation (a.r.i.), leaf water status recov- 
ered before gas-exchange variables did. At just 8h a.r.i., 
once REW was replenished (Fig. la), f experienced a large 
recovery up to values around -2.0 MPa, with higher values 
in SP than in S plants (Fig. 2a). On day 2 a.r.i., RWC recov- 
ered fully in S and SP plants, while W did on days 3 and 5 
a.r.i. in the case of SP and S plants, respectively. However, 
gs, gm, and ETR had not reached C plant values in S and 




0 2 4 7 9 11 1 2345 
Days a.w.i. Days a.r.i. 



Fig. 2. Leaf water potential (<f) (a) and RWC (b), both at 1 1 .00 GMT, 
in olive plants throughout the experiment. Values are means (±SE) 
of six replicates per treatment and date. Treatments, letters, and the 
discontinuous line as described in Fig. 1 . 
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Fig. 3. Time course during the drought period and the subsequent 
resumed irrigation of stomatal conductance to H 2 0 (g s ) (a), mesophyll 
conductance to C0 2 lg m ) (b), and net rate of C0 2 assimilation (A u ) 
(c), all measured at 1 1 .00 GMT. Values are means (±SE) of four to six 
replicates per treatment and date. Unshaded areas indicate g m data 
with a dC c /d/4 N of between 10 and 50, which is reliable according to 
Harley era/. (1992). Treatments, letters, and the discontinuous line are 
as described in Fig. 1 . 



SP plants by day 10 a.r.i. In view of these results, another 
experimental day was programmed on day 47 a.r.i. (Figs 3 
and 4), which finally confirmed the recovery in g s . In all vari- 
ables, SP plants showed a larger capacity for recovery than 
S plants insofar as, on day 4 a.r.i., g s in S and SP plants 
had recovered by 6 and 35%, g m had recovered by 35 and 
53%, and had recovered by 20 and 53%, respectively. 
Oe PIP 1.1 expression did not change significantly after the 
recovery irrigation, except for a peak observed on day 2 a.r.i. 
in S plants, similar to that found during the first stage of 
the stress period. After this, all three treatments behaved 
similarly. In contrast, OePlPl. 1 expression presented a clear 
enhanced expression after irrigation, the response being 
stronger in S than in SP plants. CA expression recovered on 
day 3 a.r.i. (Fig. 5c). 
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Fig. 4. ETR at 1 1 .00 GMT throughout the experiment. Each value is the 
mean (±SE) of six replicates per treatment and date. Treatments, letters 
and discontinuous line as described in Fig. 1 . 



Photosynthesis limitations 

During the drought period, the diffusional limitations 
(D L -S L +MC L ) prevailed over B L (Fig. 6). Two main stages 
in the time course of photosynthesis limitations could be dif- 
ferentiated in both S and SP treatments. In the initial one, 
during the first 2-4 d the drop in g s was accompanied by a 
predominant role of S L accounting for up to 50% of the T L 
in S and 35% in SP plants. In the second one, during the next 
9 d without irrigation where B L , and especially MC L , took 
the dominant role, S L lost importance. During this stage of 
severe water stress, B L was around 20% of T L . At the end of 
the drought period, all limitations in S plants were imposed 
by MC L and B L , accounting for 35 and 65%, respectively. 
Similar trends were observed in SP plants, with a minimal 
role of Si, (17%). 

During the recovery period, D h also prevailed over B h . The 
dynamics of the limitations during this period were different 
between the stressed treatments. In S plants, on 1 d a.r.i. nearly 
100% of T L was S L , and it then accounted for about 40% until 
the end of the experiment. The main difference observed 
between S and SP plants was that, in the latter, S L was not so 
exclusive on the first day of recovery, with MC L and B L play- 
ing a significant role. On the following days, T L was between 
20 and 30% lower in SP than in S plants, with a significantly 
lower MC L . All limitations had disappeared by d 47 a.r.i. 

Relationships between variables and causal models 

Normalizing the gene expression, g s , and g m by the first date 
(26 July) allowed us to eliminate the factor of different absolute 
values and concentrate on the evolution along the experiment. 
There were linear positive relationships of g m (P=0.0001) and 
g s (P=0.0001) with the expression of OePIP2.1 (Fig. 7b, e) 
and a hyperbolic one between g m and the expression of CA 
(P<0.0001) (Fig. 7c). Although OePIPl.l did not show a sig- 
nificant relationship with g m and g s , this was largely due to the 
inclusion of the peak points measured just after the imposi- 
tion of water stress or just after irrigation was resumed. If we 
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Fig. 5. Time course of the relative gene expression of aquaporins 
OePIPI .1 (a) and OePIP2.1 (b), and CA (c) in leaves of O. europaea 
throughout the experiment. Values are means (±SE) of three replicates per 
treatment and date. Treatments, letters, and the discontinuous line are as 
described in Fig. 1 . 
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Fig. 6. Quantitative limitation analysis of photosynthesis during 
the drought period and the subsequent resumed irrigation in S (a) 
and SP (b) plants. The shaded areas represent the percentage of 
stomatal (SJ, mesophyll conductance (MCJ, and biochemical (SJ 
limitation based on C plants on each date. Values are means of four 
to six replicates per treatment and date. S and SP plants were last 
irrigated on day 0 a.w.i., with daily irrigation being applied again 13 d 
later (discontinuous line). 



determined mainly by S L +MC L (Grassi and Magnani, 2005; 
Diaz-Espejo et al. , 2007); (ii) AQPs can affect H 2 0 fluxes (g s ) 
(Secchi et al., 2007ft; Laur and Hacke, 2013); (iii) AQPs can 
affect C0 2 fluxes (g m ) (Heckwolf et al., 2011; Kawase et al, 
2013); (iv) PIP1.1 and PIP2.1 interact to modify the traffic to 
the membranes and to build a tetramer that confers the basis 
of the C0 2 transport function (Zelazny et al, 2007; Otto 
et al, 2010); and (v) CA can affect g m in sclerophyll plants 
(Gillon and Yakir, 2000). Fig. 8 shows three of the most rep- 
resentative models tested. These three were consistent with 
hypotheses 1 and 2. Model 1 included the role of OePIPl.l 
in g m , reported for other PIPls in other species, but, even 
though CA was included affecting g m directly, this model 
was rejected on the grounds that P was <0.05 (i.e. it had lit- 
tle prediction capacity). Nevertheless, when in model 2, we 
changed the role of OePIPI. 1 in g m for that of OePIPI. 1 and 
included the interaction between OePIPl.l and OePIP2.1, 
simulating the tetramer structure, we obtained a satisfac- 
tory explanation for a high proportion of the variance found 
(T^O.83). However, the best fit to the data was reached with 
model 3, which completed model 2 by incorporating the role 
of CA, as a function of ^4 N , in g m (P-0.98). The inclusion of 
W (as a surrogate of plant water stress) affecting OePIPl.l 
and OePIPI. 1 reduced the proportion of the variance found 
(X 2 = 12.92; P-0.93), although the resulting model continued 
to be useful. 



discarded these values (on the bottom right corner of Fig. 7a, 
b), we obtained a significant correlation between variables in 
both cases (P=0.003, r 2 =0.43, and i><0.0001, r 2 =0.64, for g m 
and g s , respectively). 

Several alternative models were performed and studied 
with a path analysis to reveal the causal relationships that 
linked gene expression with g m , g s , and A^. Of all the pos- 
sible combinations among variables, only those that were 
physiologically sound were tested. We stated a priori the 
relationships among variables with a strong mechanistic or 
well-established and accepted empirical basis only (Shipley, 
2000). The main underlying hypotheses were: (i) is 



Discussion 

Physiological response to water stress and recovery 

Photosynthetic variables were the first to respond to irriga- 
tion withholding, being lower in stressed treatments than 
in C plants on d 2 a.w.i. (Figs 3 and 4). Leaf water status 
responded 2 d later (Fig. 2). The initial reduction in g s , larger 
in S than in SP plants, was enough to keep W at constant 
values, demonstrating the fine control of plant water status 
in olive by stomata at the initial stages of water stress (Figs 
2a, and 3a) (Fernandez et al. , 1997; Torres-Ruiz et al. , 2013a). 
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Fig. 7. Relationships between g m and gene expression (a, OePIPI .1; b, OePIP2A ; c, CA) and between g s and gene expression (d, OePIP1.1; e, OePIP2.1), 
all normalized to the first day of experiment (26 July). Values are means on each date and treatment (g m and g s : n=6; gene expression: n=3). Treatments are 
described in Fig. 1 . 



Stomatal closure has been interpreted as a plant mechanism 
to maintain xylem W well above a critical value avoiding cata- 
strophic hydraulic failure (Tyree & Zimmermann, 2002). On 
d 4 a.w.i., W was around -3.0 MPa, corresponding to a xylem 
W around -5 MPa, which is not enough to provoke impor- 
tant losses of hydraulic conductivity in olive (around 1 5-30% 
loss of hydraulic conductivity), as some works have reported 
(Ennajeh et al, 2008; Diaz-Espejo et al, 2012). However, at 
this stage, g s had nearly reached its minimum value, and there 
was no further capacity to control W, which later dropped 
to very low values of -6 MPa in the case of S treatment. The 
progress of stress in the experiment was quicker than we 
expected initially. Although we used large pots (50 1), the leaf 
area and the high atmospheric demand in our experimental 
site meant that plants consumed the available water in a few 
days. The reduction in leaf area by pruning was not enough to 
make large differences between both stressed treatments, but 
some significant differences were still found between them as 
shown in Figs 3 and 4, 2 d a.w.i. and especially a.r.i. The time 
of the day at which W and leaf gas exchange measurements 
were made were chosen carefully to try to determine the most 



significant moments: minimum W and maximum g s . The lack 
of difference in W between treatments on many days, as a 
consequence of its isohydric behaviour, did not affect the dif- 
ficulty of unravelling causal relationships, although certainly 
it would be desirable to repeat the experiment under a more 
gradual drought-stress imposition in the future. 

After rewatering, leaf water status recovered earlier than 
photosynthetic variables from drought, matching the C plant 
values on d 5 a.r.i., while ^4 N and ETR did so on d 10 a.r.i. (Figs 
2-4). Olive has a high capacity for recovery from water stress 
after rewatering, as previous studies have shown (Moriana 
et al, 2007; Torres-Ruiz et al, 20136). A similar pattern has 
been reported in tobacco (Galle et al, 2009). The value of g s 
showed the slowest recovery of all the variables studied. The 
reasons for this commonly slow recovery of g s after restoration 
of leaf water status is not yet clear, although it has been related 
to both hydraulic limitations (Brodribb and Cochard, 2009) 
and chemical limitations (Lovisolo et al, 2008). The above- 
mentioned low values of W reached at the end of the drought 
period (approx. -6 MPa) suggested that a large loss of stem 
hydraulic conductivity could be a hydraulic factor conditioning 
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Fig. 8. Examples of the alternative causal models tested using the 
d-separation method to describe causal relationships among photosynthetic 
(g s , g m , and A N ) and molecular variables (genetic expression of OePIPl.l , 
OePIP2.T , and CA). Arrows and numbers close to them indicate direct 
relationships and correlation coefficients, respectively, between variables (in 
brackets are the non-significant correlation coefficients). 



(Torres-Ruiz et al. , 20136). Other components of the hydraulic 
system could have not recovered after resuming irrigation and 
would have been limiting, like leaves or roots. However, this 
does not look to be the case, as reported recently by Torres- 
Ruiz et al. (20136). Concerning chemical limitations, several 
authors have pointed out the accumulation of abscisic acid 
(ABA) during stress acclimation as a circumstance that might 
prevent g s from fully recovering once water is available again 
and plant water status has been restored (Davies and Zhang, 
1991; Lovisolo et al, 2008). Lovisolo et al. (2008) explained 
this role played by ABA in the recovery stage as a way of con- 
trolling transpiration rate during the embolism-repairing time. 
However, the conclusions of that work, carried out in grape- 
vine, are not necessarily applicable to olive, especially as the 
former has been reported to have a great refilling capacity of 
embolized vessels (Brodersen et al., 2010). Moreover, even in 
grapevine, contradictory results have been obtained, with no 
recovery of g s once ABA accumulation was fully reversed after 
a few days (Pou et al., 2008). Finally, several authors have 
reported that recovery from water stress depends on the level 
and velocity of the stress imposition (Niinemets et al, 2009; 
Galle et al., 201 1). In this sense, our experiment showed a faster 
recovery of all physiological variables in SP than in S plants. 

Our results confirm the hypothesis that D L prevails under 
most water-stress situations in numerous species (Flexas et al. , 
2002; Keenan et al., 2010) and that B L usually appears only 
under severe water stress (Flexas et al. , 2006a; Galmes et al. , 
2007a) (Fig. 6). Although during the first 2-4 d of drought S h 
was the only prevailing limitation, as drought intensified to a 
severe level both MC L and B L increased. MC L in particular 
accounted for up to 60% of the T L , reflecting the important 
and determinant role of g m in the photosynthesis rates. This 
agrees with previous reports in similar experiments (Galmes 
et al., 2007a; Galle et al, 2011) where, under the most severe 
conditions of water stress MC L proved to be the maximum 
limitation. Once water was available again after irrigation, the 
most responsive variable determining the recovery of ^4 N was 
g m . MC L was reduced to a residual 1 5-20% in the S treatment, 
while it was nearly eliminated in the SP treatment (Fig. 6). 
The faster recovery of photosynthetic variables in SP com- 
pared to S plants translated into lower values of T L in SP 
than in S plants. A good correlation between g s and g m was 
found for most of the studied period with a g m :g s ratio (both 
on C0 2 basis) of 2.1. Considerable evidence of the close rela- 
tionship between g s and g m has been reported (Flexas et al. , 
2008; Warren, 20086), but this relationship has been regarded 
as the result of a tight covariance of these two resistances 
in the C0 2 pathway to the leaf rather than as a cause-effect 
relationship. In fact, several works have indicated that the 
relationship between g s and g m can be modified according to 
changing environmental conditions (Flexas et al, 2009; Galle 
et al, 2009 Perez-Martin et al, 2009; Flexas et al, 2013a). 



the slow and partial recovery of g s when water was applied 
during the recovery stage. However, even at these values of W, 
olive is not expected to lose much of its hydraulic conductiv- 
ity: between 35 and 40% according to Ennajeh et al. (2008), 
explaining the ability of this species to tolerate severe drought 



Gene expression ofAQPs and CA during 
drought and recovery 

The expression patterns during stress treatments were differ- 
ent between the two AQPs studied and were more irregular 
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than for CA (Fig. 5). This impression was shown by the peak 
responses measured at the beginning of both stress and recov- 
ery periods in OePIPl.l. The peaks at the onset of drought 
have been reported in similar studies (Yamada et al, 1997; 
Galmes et al, 2001b; Pou et al, 2013), where the upregula- 
tion of AQP expression, including PIP1 and PIP1, was inter- 
preted as a mechanism to promote water movement inside 
leaves via symplast by increasing membrane permeability to 
water when this is less available for the plant. The peak of 
expression in OePIPl.l at the start of the re-irrigation period 
could explain the fast recovery of leaf water status by the pos- 
sible role of AQPs in the xylem refilling of parenchyma cells 
(Kaldenhoff et al , 2008). This peak would occur in a moment 
of high water availability in soil, which would reduce risks of 
massive losses of water from leaves (Secchi et al, 2007a). On 
the other hand, the downregulation in OePIP2.1 expression 
at the end of drought period may encourage cellular water 
conservation during periods of water stress by reducing mem- 
brane water permeability and limiting loss of cellular water 
(Secchi et al, 2007a). Probably, in order to maintain a suit- 
able water status under abiotic stress, both increased water 
transport via AQPs in some tissues and reduced in others are 
required (Secchi et al, 2001b) because one of the main role 
of AQPs is to maintain homeostasis and water balance under 
water-stress conditions (Tyerman et al, 2002). 

In addition to their role in water transport, there is pub- 
lished evidence of the putative role of AQPs in C0 2 trans- 
port and g m (Uehlein et al, 2003; Kawase et al, 2013). 
Accordingly, we found a linear relationship between 
OePIPl.l expression and g m , but no relationship between 
OePIPl.l expression and g m (Fig. 7). Although it is found 
that most members of PIP2 act as water channels (Tsuchihira 
et al, 2010) while PIP1 members are involved in a C0 2 
transport function (Flexas et al, 2006ft; Heckwolf et al, 
2011), Hanba et al. (2004) found differences in g m induced 
by a PIP2. On the other hand, it has been reported that PIP1 
and PIP2 AQPs form heterotetramers, which modify their 
function as H 2 0 or C0 2 membrane transport facilitators 
depending on their composition (Otto et al, 2010); hence, 
both PIP1 and PIP2 may be inextricably linked with each 
other and with g m . On the other hand, OePIP2.1 showed a 
similar degree of correlation with g s , as already observed for 
PIP2.1 of grapevines (Pou et al, 2013). In summary, while 
it is unclear from the data whether the expression of specific 
AQPs is related to regulation of both g s and g m , or whether 
they just operate by affecting one of the two conductances 
and a tight co-regulation between g s and g m provokes an 
apparent correlation of the other with PIP expression, the 
results strongly suggest that AQP expression is involved in 
setting diffusional limitations to photosynthesis in olives 
under water stress and recovery. This same link provokes an 
additional tight correlation between the two apparently less 
related parameters of g m and the leaf hydraulic conductance 
(Flexas et al, 2013ft). 

The CA expression patterns in stressed treatments through- 
out the experiment were smoother than AQP patterns and 
were similar in S and SP plants (Fig. 5c), closely tracking those 
observed for ETR (Fig. 4). The role of this stroma CA in g m 



regulation was supported by the relationship found between 
CA expression and g m (Fig. 7c). The hyperbolic shape of this 
relationship agrees with the consideration of CA as not being 
limiting to photosynthesis due to the large amount of CA pre- 
sent in leaves, although this can be species dependent (Makino 
et al, 1992). Although the role of CA in regulating g m is more 
controversial than the evidence found for AQPs (Flexas et al , 
2012), some authors have reported that this role can be spe- 
cies dependent, gaining importance when g m is low, as hap- 
pens in sclerophyll species (Gillon and Yakir, 2000). Therefore, 
this would also be the case for olive, a sclerophyllous species 
with a thick cell wall (Bacelar et al. , 2004; Marchi et al. , 2008). 
Following the results of Gillon and Yakir (2000), olive would 
have to counterbalance the low conductance at the cell-wall 
sites with a larger conductance at the chloroplast sites, which 
undoubtedly would make g m more dependent on CA. 

To analyse the causal nature and the structure of the 
relationships found here (Fig. 7b, c), we considered several 
causal models with physiological sense that fitted the data. 
The most robust structure with the highest significant levels 
in many of the combinations tested was that obtained in 
model 3 (Fig. 8). This model showed that there could be a 
direct role of OePIPl.l and CA expression in g m regula- 
tion, as well as an indirect role of OePIPl.l. The interac- 
tion between both AQPs is justified in accordance with the 
conclusions by Otto et al. (2010), who suggested that PIP1 
and PIP2 interact to form a heterotetramer. Additionally, 
the model considers that both PIPs are directly involved 
in the regulation of g s , that both S L and MC L determine 
yl N , and that CA is in turn regulated by A N . OePIPl.l and 
OePIPl.l in olive have been described to regulate mem- 
brane water permeability in leaves (Secchi et al. , 2007ft), a 
fact that explains their effect on g s . CA was tested in an 
independent model (not shown) to regulate g m directly with- 
out the role of AQPs, this model being able to explain only 
a small amount of the data variance (x 2 =22.53; degrees of 
freedom=16; P-0A26). This supports the hypothesis that 
the role played by CA in the regulation of g m is not central. 
In any case, in addition to the implications in sclerophyl- 
lus leaves mentioned above (Gillon and Yakir, 2000), the 
inclusion of CA as a function of ^4 N is justified because the 
regulation of CA gene expression is related to the leaf pho- 
tosynthetic activity by the C0 2 inside leaves (Hoang and 
Chapman, 2002; Fabre et al. , 2007). CA has been reported 
to be closely related to Ribulose-l,5-bisphosphate car- 
boxylase/oxygenase (Rubisco) activity, and therefore to A N 
(Makino et al, 1992; Price et al, 1994). Furthermore, the 
time course of CA expression and A N were quite similar 
throughout the experiment (Figs 5c and 3c) and recent stud- 
ies have demonstrated the role of CA as an upstream regu- 
lator of C0 2 -controlled stomatal movements in guard cells 
(Hu et al, 2010). While the validation of a causal model 
does not imply that the hypothesis is true, it does imply that 
the hypothesis is plausible given the empirical data (Shipley, 
2000). C; has long been proposed as a potential signal used 
by the plant to maintain the equilibrium among g s , g m , and 
A N (Ethier et al, 2006; Flexas et al, 2009). Under water 
stress, downregulation of g s , g m , F cmax , and ETR could 
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occur to keep Q under relatively constant values during the 
acclimation process. If this hypothesis is correct, CA would 
be in an optimal position for regulation of g m , as it uses 
as Q as a substrate and influences its diffusion inside the 
cell (Terashima et al. , 201 1). Although gene expression may 
not necessarily reflect protein function because of the post- 
transcriptional regulation (Maurel, 2007; Heinen et al, 
2009; Laur and Hacke, 2013), the good relations of g m with 
expression of OePIP2. 1 and C A and the good fitting of the 
causal model proposed allow us to go beyond already pub- 
lished literature and to infer a direct role of AQPs and CA 
in g m regulation, although certainly the mechanistic basis 
for such a role remains to be elucidated. 

Conclusions 

The main objective of this study was to provide a step for- 
ward in the putative role played by AQPs and CA in the 
regulation of g m under water stress. Despite the use of the 
path analysis to infer some potential causal effects, a novel 
aspect of this work is the execution of the experiment out- 
doors under natural conditions. Photosynthetic variables 
in olive trees responded faster than water status under 
drought but recovered much later. Despite the severe stress, 
D L prevailed over B h throughout the experiment. The role 
of MC L varied throughout the experiment, prevailing over 
the rest of the limitation components as stress intensified 
and being present even during the recovery phase. This 
study supports the hypothesis that the regulation of g m is 
regulated mainly by AQPs and that both OePIPl.l and 
OePIP2.1 are likely to interact to exert a significant effect 
on g m . The data shown in this study reveals novel evidence 
on the putative role of CA in the regulation of mesophyll 
conductance to C0 2 , which, although small if compared 
with that of AQPs, is justified by the sclerophyllous nature 
of olive leaves. Due to the tight correlation between AQPs 
and both g s and g m , evidence of co-regulation between the 
two conductances and the fact that H 2 0 and C0 2 fluxes in 
leaves share a part of their pathways, facilitated by AQPs, 
it is very difficult to relate genes with functions in an une- 
quivocal manner. As far as we know, there are few than 
expected studies showing this sort of correlation and, of 
course, no reverse genetics approach is possible to date with 
woody species like olive trees. Hence we believe that our 
approach, although certainly limited, is useful for a bet- 
ter understanding of plant physiological response to water 
stress, and is perhaps the only one currently available for 
olives. New biochemical probes and techniques are there- 
fore required in order to quantify the activity of AQPs and 
their genetic expression changes, as well as to follow their 
spatial distribution in cells. 

Supplementary data 

Supplementary data are available at JXB online. 

Supplementary Fig. SI. Picture showing the aspect of the 
experimental set-up: 5-year-old trees in 50 1 pots. 



Supplementary Fig. S2. Relationship of the quantum yield 
of photosystem II (OPSII) and OCO 2 [(^ N +.R d )/PPFD]. 

Supplementary Fig. S3. Relationship between V cmax and 
ETR used in the quantitative analysis of photosynthesis limi- 
tations throughout the experiment. 

Supplementary Fig. S4. Maximum-likelihood tree of 
Oka europaea carbonic anhydrases (CA) with Arabidopsis 
thaliana CA. 

Supplementary method. Assessment of limitations to 
photosynthesis. 
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